Effect of pressure on the phase behavior and structure of water confined between 
nanoscale hydrophobic and hydrophilic plates 
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We perform systematic molecular dynamics simulations of water confined between two nanoscale 
plates at T = 300 K. We investigate the effect of pressure (-0.15 GPa < P < 0.2 GPa) and 
plate separation (0.4 nm < d < 1.6 nm) on the phase behavior of water when the plates are 
either hydrophobic or hydrophilic. When water is confined between hydrophobic plates, capillary 
evaporation occurs between the plates at low enough P. The threshold value of d at which this 
transition occurs decreases with P (e.g., 1.6 nm at P « —0.05 GPa; 0.5 nm at P 0.1 GPa), 
until, at high P, no capillary evaporation occurs. For d ~ 0.6 nm and P > 0.1 GPa, the system 
crystallizes into a bilayer ice. A P — d phase diagram showing the vapor, liquid, and bilayer ice 
phases is proposed. When water is confined by hydrophilic (hydroxylated silica) plates, it remains 
in the liquid phase at all P and d studied. Interestingly, we observe, for this case, that even at 
the P at which bulk water cavitates, the confined water remains in the liquid state. We also study 
systematically the state of hydration at different P for both kinds of plates. For the range of 
conditions studied here, we find that in the presence of hydrophobic plates the effect of P is to 
enhance water structure and to push water molecules toward the plates. The average orientation of 
water molecules next to the hydrophobic plates does not change upon pressurization. In contrast, in 
the presence of hydrophilic plates, water structure is insensitive to P. Hence, our results suggest that 
upon pressurization, hydrophobic plates behave as 'soft' surfaces (in the sense of accommodating 
pressure-dependent changes in water structure) while hydrophilic walls behave as 'hard' surfaces. 
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I. INTRODUCTION 

Confining geometries that contain thin water films are 
common in biology, geology and engineering. Examples 
include ion channels lL ,1 H S Sll, mesoscopic sur- 
factant assemblies 0, , mineral inclusions , zeolites 
[lol |. and microfluidic technologies |^. Understanding 
the changes in water structure, dynamics and thermo- 
dynamics due to interactions with confining surfaces is 
therefore of interest in a wide variety of scientific prob- 
lems and technical applications, such as the inhibition of 
corrosion T^l, heterogeneous catalysis, the design of su- 
perhydrophobic surfaces ^3 ' the ascent of sap in plants 
[3, the function of biological membranes and the 
development of 'lab on a chip' technologies p^Sj . 

The behavior of water near hydrophobic surfaces has 
attracted considerable attention, both on account of the 
intrinsic scientific interest of the thermodynamic prob- 
lem (e.g., |ililLa8(,a9t^i21»^2^2415, 26J), as 
well as its relevance to water permeation through mem- 
brane channels [^11,1^1^ sl I27II and carbon nanotubes 



[28, 29] , micelle formation 
ity of membrane proteins 




and the amphiphatic- 
33 . Hydrophobic in- 



teractions are also important in phase transfer catalysis 
Issl l. chemical self-assembly of macroscopic objects js^ . 



and capillary evaporation [s^ |3^ [s^ ■ 
More than thirty years ago, Stillinger ^ 



argued that 



when a non-polar solute is sufficiently large, the water- 
solute interface resembles that between water and its 
own vapor. Lum et al. proposed a theoretical ap- 
proach to describe the crossover between the solvation 
of small and large hydrophobic objects. Bulk thermo- 
dynamics and entropic contributions to the free energy 
of solvation are dominant in the former case. Interfa- 
cial thermodynamics and enthalpic contributions to sol- 
vation free energy play a key role in the latter case. 
The crossover length scale between these two regimes, 
w 1 nm at ambient conditions, is a sensitive function 
of pressure and the presence of co-solutes [s^. From 
a microscopic point of view, this length scale is related 
to the different manner in which water molecules ar- 
range around solutes of different sizes. In the presence of 
small hydrophobic solutes, water molecules rearrange in 
such a way that the solutes can be accommodated into 
water's hydrogen-bond (HB) network. Thus, small hy- 
drophobic solutes are caged in clathrate-like structures 
m E EJl with HB vectors (defined as the four tetra- 
hedral vectors pointing outward from the oxygen atom 
of each water molecule along the two oxygen-hydrogen 
bond donor, and the two "lone pair electrons" acceptor 
directions) avoiding pointing toward the non-polar so- 
lute. As the solute size increases (to « 1 nm at ambient 
conditions) ji^, 113, 0, j such clathrate-like structure 
cannot be maintained: HBs break, inducing a restruc- 
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turing of the HB network. In this case, the hydration 
structure shows an orientational inversion with respect 
to the clathrate-hke structure, with HBs pointing toward 
the solute IH El H gllil . 

When two hydrophobic surfaces approach each other, 
a confined water film can be thermodynamically destabi- 
lized with respect to the vapor at small enough separa- 
tions ,21, • A plausible interpretation of the 
experimentally observed long-ranged attractions between 
hydrophobic surfaces (e.g., 10-100 nm) |H1 El i| 
involves dewetting of the inter-plate region. As shown 
by Berard et al. [5J| , such dewetting is not particular to 
water, but arises as a general consequence of weak fluid- 
wall interactions relative to fluid- fluid interactions (2^. 
The generality of confinement-induced drying notwith- 
standing, interest in this phenomenon has focused pre- 
dominantly on the case in which the confined fluid is 
water and the confining surfaces are hydrophobic. Sev- 
eral coinputational studies of this situation exist (e.g., 

111 H iHi ii il il m m in ) . 

Kauzmann's influential review '62^ first pointed out the 
importance of hydrophobic interactions in protein fold- 
ing. It is generally accepted that the formation of a core 
of hydrophobic residues from which water is partially 
excluded is a dominant force in the folding of g lobular 
proteins [llillMlillillillllillTaEll How- 
ever, the surface of a protein is a geometrically, chemi- 
cally, and electrically heterogeneous object, and the dis- 
tribution, structure, and dynamics of water near such 
a surface is expected to be quite different from that 
found near idealized surfaces, such as perfectly hydropho- 
bic walls. In fact, computer simulations have clearly 
shown a range of water behavior in response to surface 
heterogeneity. Cheng and Rossky 72] studied the hy- 
dration of the polypeptide melittin. They found that 
clathrate-like structures dominate near convex surface 
patches, while near flat regions the hydration shell fluc- 
tuates between clathrate-like and less ordered (inverted) 
structures. The collapse of a two-domain BphC enzyme 
showed no dewetting in the inter-domain region when 
the units were brought together down to a distance of 
0.4 nm j?!. A dewetting transition was observed only 
when electrostatic interactions were turned off Ws^. In 
contrast, dewetting has been observed in the collapse of 
of the melittin tetramer [t^; however, even single muta- 
tions were shown to prevent drying. 

The above protein folding examples illustrate the im- 
portance of understanding water behavior near complex, 
heterogeneous surfaces. Materials science and engineer- 
ing applications involving chemically nanopatterned sub- 
strates (e.g., [t^) likewise will require an understanding 
of water structure near heterogeneous interfaces. Molec- 
ular simulations should prove extremely powerful in ad- 
dressing this problem, because they allow exquisitely sen- 
sitive control of surface geometry and chemistry. Ques- 
tions such as the manner in which the characteristic 
length of hydrophobic patches and the pattern of their 
distribution on a hydrophilic surface affect confined wa- 



ter structure and dynamics, or the effect of surface chem- 
istry in promoting or preventing bulk cavitation follow- 
ing dewetting of a confined region are ideally suited for 
computational scrutiny. The recent work of Koishi et 
al. (59j is an excellent example of the valuable insights 
that can be obtained by molecular simulations regarding 
water behavior near complex surfaces. 

This paper is the first report of an ongoing compu- 
tational investigation aimed at understanding the influ- 
ence of surface heterogeneity and bulk thermodynamic 
conditions (pressure, temperature) on the structure, dy- 
namics, and thermodynamics of confined water. Here we 
establish the 'base case' conditions by comparing water 
structure and phase behavior at 300 K when confined 
by purely hydrophobic and purely hydrophilic nanoscale 
surfaces (0.4 < d < 1.6 nm). Li et al. |53| have recently 
studied hydration and dewetting near hydrophilic and 
hydrophobic plates (but with different chemistry than 
the ones we consider here), at a single pressure. An 
important aspect of the present work is that we per- 
form such a comparison over a broad range of pressures 
(—0.15 < P < 0.2 GPa). In subsequent papers we will 
report results on water confined between heterogeneous 
plates with engineered 'patchiness'. 

Our computational study should be placed in the con- 
text of a considerable body of experimental work aimed 
at studying systematically the properties of water con- 
fined by hydrophobic and hydrophilic surfaces. The 
forces between hydrophobic surfaces in aqueous solution 
have been measured with the surface force apparatus 
(e.g., m m m m 0). such studies revealed the 
existence of long-range attractions between hydrophobic 
surfaces, over distances of the order of 100 nm, and in- 
vestigated the origin of this phenomenon. Bridging cavi- 
ties and microbubbles were linked in several I sol 153 . l53j| , 
but not all 51] of these studies, to the long-ranged at- 
traction. Zhang et al. [t^ investigated the response to 
shear stresses of water confined between adjoining sur- 
faces, one hydrophobic and the other hydrophilic (Janus 
surface) . They found that the competing effects of these 
surfaces gave rise to unusually noisy responses to shear. 
Jensen et al. [SI] used X-ray reflectivity to study the 
contact region between water and an extended paraffin 
surface. They found that drying was confined to a very 
narrow (w 1.5 nm) region. Ruan et al. [t^ investigated 
the structure and dynamics of water near a hydrophilic 
surface using ultrafast electron crystallography. They ob- 
served coexistence between ordered surface water (up to 
1 nm thick) and crystalline islands. Reference to particu- 
lar aspects of the studies summarized above will be made 
throughout this paper, stressing, whenever possible, the 
relationship between experimental observations and our 
calculations. 

This paper is organized as follows. In the next section 
we describe the simulation details. The results for hy- 
drophobic plates are presented in Sec. IIIII calculations 
on water confined by hydrophilic plates are reported in 
Sec. IIVI We summarize our main conclusions in Sec. IVl 
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II. SIMULATION DETAILS 

We perform molecular dynamics (MD) simulations in 
the NPT ensemble. The system is composed of = 3375 
water molecules in which are inmersed two identical finite 
three dimensional nanoscale plates. The temperature is 
fixed at T = 300 K by using a Berendsen thermostat 
[t^. The pressure ranges from P = —0.15 GPa to P = 
0.2 GPa, and is controlled by coupling the system volume 
to an external bath at P (analogous to the Berendsen 
thermostat [t^I. 

We simulate a cubic sample using periodic boundary 
conditions along the three directions. Water molecules 
are modelled with the SPC/E pair potential js^- The 
plates are introduced symmetrically about the center of 
the box such that they are parallel to the x — y plane and 
equidistant from the z = plane. The plate dimensions 
(3.215 X 3.217 X 0.866 nm^) are smaller than the box size 
(which varies with P but always exceeds 4.85 nm under 
the conditions investigated here). Figures ^a) and (b) 
show a fully hydroxylated silica plate (one of the two 
kinds of plates used in this study) and Fig. QJc) is a 
cross section, showing the water molecules and the two 
cavities where the plates are located. The plates are fixed 
in space throughout. 

We consider two kinds of plates: hydrophobic or hy- 
drophilic. Their common underlying structure corre- 
sponds to four layers of Si02 [see Fig. ^ ^b)! reproducing 
the (1.1.1) octahedral face of cristobalite |8lLl82 | [see Fig. 
n^a)]. The unit cell of Si02 is idealized as a perfect tetra- 
hedron with O — O and Si — O distances of 0.247 nm and 
0.151 nm, respectively. The hydrophobic plates consist of 
the above-described structure; however, the 'S'i' and 'O' 
atoms interact with water molecule O atoms exclusively 
via a Lennard- Jones potential (e and a parameters are 
given in Table P) . Each hydrophobic plate is composed 
of 674 atoms. 

The hydrophilic plates correspond to fully hydroxy- 
lated silica and are obtained by attaching a hydrogen 
atom to each surface oxygen atom on the four plate sur- 
faces [see Figs, ^d) and (e)]. The O — H distance is 
chosen to be the same as in the SPC/E model, i.e. 0.1 
nm. The Si and O atoms are located in fixed positions 
(as in the hydrophobic plates) but the H atoms on the 
surface are able to move with fixed bond lengths and 
bond angles; each H atom can reorient in a circle. Such 
circular motion occurs in a plane parallel to the plate at 
a distance 0.033 nm away from the O atom plane of the 
Si — O — H groups. The resulting Si — O — H angle is 
109°. 27 [see Fig. ^e)]. Each hydrophilic plate is com- 
posed of 778 atoms. Only the polar Si~0 — H groups at 
the surface carry partial electric charge [charge values are 
given in Table H] . Electrostatic interactions are treated 
using the Ewald sum method with a cutoff distance of 
0.79 nm and parameters rrimax = 5'^ (for the number of 
vectors in the reciprocal-space sum) and a = 0.4 (for 
the width of the screening-charge Gaussian distribution) 

Ulii. 



The vibrations of the Si and O atoms have not been 
taken into account in this work. A computational study 
of water in single-wall carbon nanotubes indicates that 
flexibility can affect a channel's apparent hydrophobicity 
[29] . On the other hand, a molecular dynamics study 
of water droplets on graphite revealed a negligible effect 
of substrate vibrations on the contact angle [s^. The 
lattice constant for the cristobalite structure used in this 
work is 0.494 nm. Invoking a Lindemann-type estimate, 
the amplitude of individual atomic vibrations should be 
considerably smaller than the melting threshold (10% of 
the lattice constant), i.e. 0.05 nm in our case. This upper 
bound on the magnitude of substrate atomic vibrations 
is an order of magnitude smaller than the closest inter- 
plate separation considered in our work. Thus, while the 
effect of crystal vibrations on the structure and dynamics 
of confined water deserves careful attention and will be 
explored in future studies, we believe that the present 
rigid wall base case is a reasonable starting point for our 
investigations. 

The pressure is calculated using the virial expression, 
taking into account the fact that some atoms in the plates 
are fixed and others can move [s^ . We also use the 
link cell and neighbor list methods to calculate the pair 
interactions [s^. 

We perform simulations for different time intervals de- 
pending on the kind of plates and their separation. Simu- 
lation times are indicated in Table HTl One of the reasons 
we chose these simulation times is that when there is 
no capillary evaporation or crystallization, we find that 
quantities such us total energy and volume are constant 
for t > 50 ps. Thus, the first 50 ps of the simulation are 
discarded and the rest of the simulation is used for data 
acquisition. Moreover, the correlation time (obtained 
from the intermediate scattering function) in bulk wa- 
ter simulations using the SPC/E model at T = 300 K 
and p — 1.0 g/cm'^ is approximately 3 ps 's^. Thus, our 
simulations appear long enough to avoid effects of the 
equilibration process on our calculations. 



III. RESULTS: HYDROPHOBIC PLATES 
A. Capillary evaporation: effect of pressure 

In this section, we describe the effects of pressure P and 
plate-plate separation d on the phase behavior of water 
confined between hydrophobic plates at T = 300 K. To 
consistently define d, we use the distance between the 
planes containing the hydrogen [H) atoms on the inner 
surfaces of two hydrophilic plates. In the case of hy- 
drophobic plates, where no hydrogen atoms are present, 
the same planes (where the H atoms would hypothet- 
ically be located) are used to define d. As previously 
noted, the plane containing the H atoms is located at a 
distance 0.033 nm from the plane containing the oxygen 
atoms of the silanol groups. 

Figure[21summarizes the results of our MD simulations 
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for d > 0.5 nm. Two observations are most relevant 
from Fig. [3 (i) the threshold separation distance dth re- 
quired for capillary evaporation decreases with increasing 
P, and (ii) for d « 0.6 nm a bilayer crystal is obtained 
at high P. Simulations for d < 0.4 nm at P = 0.2 GPa 
and GPa show that molecules are trapped interstitially 
between the plate atoms. The water molecules remain 
in a plane and are not able to diffuse due to the lack of 
space between plates. 

The decrease in dth with increasing P can be under- 
stood as follows. Equating the grand potential of the 
confined liquid and the confined vapor, for sufficiently 
large surfaces js^ [H^] , 

- P Adth + 2 A-f^i^- P* A dth + 2 A 7„„ (1) 

whence 



Mh 



2 A7 
{P-P*) 



(2) 



In the above equations, A is the plate surface area; 7^; 
and 7^1, are the wall-liquid and wall- vapor interfacial ten- 
sions (A7 = 7u,; —Jwv)] P is the bulk pressure; and P* is 
the equilibrium vapor pressure at the given temperature. 
It follows from Eq. (0) that dth decreases with increasing 
P, as manifest in Fig. |21 For an incompressible liquid, 
P — P* = pi {fii — fiv), and, therefore |i21j.i35|. 



Hh 



2 A7 



Pi {fJ-l - Pv) 



(3) 



where pi and p^ denote the chemical potentials of the 
liquid at P and of the vapor (and liquid) at P*. In the 
above derivation we have assumed that the plate char- 
acteristic size, « A^l"^ , is large enough so that one can 
neglect the finite lateral size of the confined region (i.e., 

dthiA^i^ < 1) m. 

A sequence of snapshots showing a typical dewetting 
process is given in Fig. |21 We observe that the density 
rapidly becomes rarified on a time scale of « 10 ps, and 
then the solvent retreats on a somewhat longer time scale 
(« 100 ps). At P = GPa and d < 0.6 nm, the vapor 
cavity is limited to the area of the plates, i.e. the water 
surrounding the plates remains in the liquid state [Fig. 
ISJ^d)]. However, at the lower pressure P = —0.05 GPa 
and rf < 1.2 nm, the bubble originally formed between 
the plates expands to the bulk liquid, inducing cavita- 
tion of the whole system. In other words, a heteroge- 
neous nucleation event occurs induced by the plates. For 
P < —0.1 GPa and d < 1.6 nm, we observe a simulta- 
neous cavitation inside and outside of the confined space 
indicating that, for the present water model, the liquid 
phase is unstable for P < —0.1 GPa. We note that when 
capillary evaporation occurs, no water molecule is ob- 
served between the plates. Thus, practically speaking, 
what we call 'vapor phase' is here a vacuum. 

Observations consistent with this picture have 
been reported in measurements of the force between 
hydrophobically-coated mica surfaces immersed in water 



[53. Specifically, vapor cavities formed spontaneously 
when mica surfaces coated with a double-chain cationic 
fluorocarbon surfactant approached to distances between 
1 and 4 nm '5^. However, cavities only formed after sep- 
aration from contact when the mica surfaces were coated 
with a double-chain cationic hydrocarbon surfactant |53j | . 
Both the separation at which cavities formed and the 
magnitude of the force between the hydrophobic surfaces 
were found to be very sensitive to the specific substance 
used to coat the mica surface [F^, E^l- The measured 
width of the dewetting layer formed at an extended paraf- 
fin surface in water (« 1.5 nm) is also consistent with the 
present calculations [55i |. 

Our results on capillary evaporation are in agreement 
with previous grand canonical ensemble Monte Carlo 
simulations of water using the SPG model [s^. That 
work shows that capillary evaporation between smooth 
planar hydrophobic walls occurs at walls separations 
D w 1.27 nm and £» « 0.9 nm at P = GPa and 
P w 0.1 GPa, respectively. These distances correspond, 
in our case, roughly to d 0.79 nm and d 0.42 nm 
(see Fig. 121. 

The presence of a crystal at high P has not been dis- 
cussed so far in the context of capillary evaporation. 
However, there is considerable evidence from MD sim- 
ulations showing that confined water may crystallize at 
high P (see e.g. 90] ). Koga et al. found crystallization in 
MD simulations of the TIP4Pjnodcl at T « 230 - 300 K 
and P = 0.5 - 1.0 GPa 1119^ ^3- The ice structure 
that we find is similar to that found in [9ll l9^ and re- 
sembles none of the structures of the existing ice phases 
in bulk water J nor those found in metallic or hydropho- 
bic pores mllEllllgll. Figure Ua), a top view of the 
bilayer ice, clearly shows a hexagonal lattice, while Fig. 
EJb) suggests that each layer is almost flat, and that the 
two layers are in registry. 

To confirm that confined water crystallizes at high P 
and d « 0.6 nm, we calculate the mean square displace- 
ment (MSD) parallel to the plates at P = 0.2 GPa for 
different values of d. Figure [S{ a) shows the MSD for a 
400 ps simulation averaging over molecules (we do not 
average over starting times). For d = 0.4 nm, we find 
that MSD K, consistent with the earlier observation 
that molecules are trapped between the plate atoms. For 
d ~ 0.5, 0.8, and 1.0 nm the MSD increases monoton- 
ically with time and, for a fixed time, with d. How- 
ever, the MSD for d = 0.6 nm shows a fast increase for 
^ ^ 50 ps and then approaches an asymptotic value of 
w 0.2 nm. Snapshots of the system during the first 50 ps 
show that the molecules in the confined space reorganize 
to form the crystal. For 50 ps < t < 300ps, the resulting 
ice shows defects that disappear with time. The trans- 
formation is almost over for t k, 300 ps. In the absence of 
crystallization, one would expect the MSD at long times 
for d = 0.6 nm to fall between that corresponding to 
d = 0.5 nm and d — 0.8 nm. Instead, the MSD curve 
for d — 0.6 nm shows a plateau, as for d = 0.4 nm, after 
about 250 ps, consistent with the view that the system 
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crystallizes. 

The crystallization of confined water is also confirmed 
by structural properties such as the radial distribution 
function (RDF) parallel to the plates, gxyir), and the 
probability density function to find a molecule located at 
z between the plates, Pn{z). FigureEIb) shows gxy{^) for 
the ice {d = 0.6 nm) and for a bilayer liquid [d = 0.8 nm). 
The oscillations in gxy{f) for d — 0.6 nm are a clear sign 
of crystallization. Such oscillations are not present in the 
liquid phase. We note that the slow decay of gxy{r) for 
large r is due to the finite size of the plates. This finite 
size effect can be removed by taking the ratio of gxy{r) 
for the crystal to that of the liquid (see inset). The per- 
sistence of oscillations is then clear. Figure |SJc) shows 
the Pn{z) for the ice (d = 0.6 nm) and for a comparable 
liquid (d = 0.8 nmV While in both cases water molecules 
form two layers [Og, the peaks of Pn{z) are much larger 
at d = 0.6 nm than at d = 0.8 nm. P„(z) was calcu- 
lated from a histogram of water molecules in 78 slices 
into which the distance between plates was divided. It is 
defined so that / Pn{z) dz = 1. 

We note that no spontaneous crystallization in bulk or 
confined water has been reported in MD simulations us- 
ing the SPC/E model. Crystallization using the TIP4P 
[99.] and TIP5P models has been obtained lOOj . In our 
case, we note that the location of water molecules belong- 
ing to each of the two ice layers is highly correlated with 
the plate structure, suggesting that the crystal we ob- 
tain is templated by the plate. Such substrate-templated 
crystallization of interfacial water has been observed on 
hydrophilic surfaces |78|| . Figure EJa) is a schematic di- 
agram showing both the hexagons formed by the water 
oxygen atoms and the atoms at the plate surface. The 
atoms at the plate surface are arranged in tetrahedra 
pointing either into or out of the page of the figure. These 
two kinds of tetrahedra alternate, forming rings com- 
posed of six tetrahedra [see also Fig. ^a)] . Two rules de- 
termine the structure of a given hexagon of the ice layers: 
(i) a water oxygen atom is located on top of each tetrahe- 
dron pointing into the page of the figure; and (ii) a water 
oxygen atom is located on top of the center of the hexag- 
onal rings formed by the tetrahedra at the plate surface. 
As a result, the ice layers are formed by perfect hexagons 
with an OO distance of 0.29 nm. Interestingly, this is the 
same distance separating the two ice layers [see Fig.|5fc)]. 
Thus, the OO distance between nearest and next near- 
est neighbors are 0.29 nm and VS x 0.29 nm = 0.49 nm, 
respectively. These values are in agreement with the lo- 
cation of the first two maxima of gxyir) [see Fig. Elb)] 
which are located at ~ 0.29 nm and « 0.51 nm, respec- 
tively. 

The bilayer ice is fully hydrogen-bonded and each 
molecule has four HBs, as in ordinary ice. HBs occur ei- 
ther between molecules in the same ice layer or between 
both layers (i.e., there are no HB pointing toward the 
plates) . The two water layers are linked by HBs where a 
molecule of one layer shares one of its H atoms with the 
nearest oxygen atom of the other layer [see Figure Etb) 



and (c)]. We note that the HOH angle in the SPC/E 
model is 109.47°, i.e., the tetrahedral angle. However, 
in the bilayer ice structure, the OOO angles (formed be- 
tween three nearest neighbors) are either 90° or 120°. 
Therefore, in general, the OO directions (between near- 
est neighbors) differ slightly from the OH directions of a 
given water molecule. 

In refs. [lOll [lo2, it was found that TIP4P water 
within carbon nanotubes can crystallize forming a tube 
of square, pentagonal, or hexagonal ice, depending on the 
nanotube radius. The structure of the carbon nanotube 
hexagonal ice seems to be the same as that found in the 
present work. Interestingly, the same requirements [that 
the confinement geometry allows: (i) next-nearest oxygen 
atoms to be loca ted roughly at the bulk 0-0 distance 
(« 0.28 nm) | 102| : and (ii) water molecules to have four 
hydrogen bonds] are necessary for ice formation in our 
case and in the carbon nanotubes. 

Simulations of water betwee n cry stalline infinite walls 
using the TIP5P model ^ Eol [1041, and between 
smooth infinite walls using the TIP5P 105j and TIP4P 
models j92i also show the formation of crystals at T = 
300 K. Both mon olayer and tri-layer ices have been re- 
ported |104 llOSj . As indicated in Fig. 13 our simula- 
tions at P = 0.2 GPa show no sign of monolayer ice at 
small d nor of crystallization to other solids at larger d 
(corresponding to n-layer crystals, n > 3). Small vari- 
ations in d (< 0.1 nm) were found to have a profound 
effect on the appe arance and disappearance of confined 
ices m Eol m^. In Fig. [3 the values of d simulated 
at a given P are separated by Ad > 0.1 nm. Thus, it is 
possible (but unlikely) that these ices can be found for 
confined SPC/E at T « 300 K when using a smaller Ad. 

Figure |21 also shows a schematic confined phase dia- 
gram superimposed on the simulation data points. The 
bilayer crystal and liquid phases are indicated, together 
with an estimate of the dewetting transition line. At 
T = 300 K, we find that the dewetting transition line ex- 
tends down to d = 0.5 nm, suggesting that it merges with 
the region (defined by d < 0.4 nm) for which molecules 
gets trapped between the plate atoms. It is interesting 
to consider whether for a higher T the dewetting tran- 
sition line might end in a critical point. The relation- 
ship between the dewetting locus and water's ordinary 
vapor-liquid critical point will be investigated in future 
studies at higher temperatures. It is also possible that 
the bilayer crystal region in Fig. |21 expands down to the 
dewetting transition line (see red dotted line in the fig- 
ure). In this case, there would be a triple point where the 
vapor phase coexists with the liquid and crystal phases. 
Our simulation at P = 0.05 GPa and d — 0.6 nm does 
not show the presence of a crystalline structure. Instead, 
we observe a competition between a disordered (liquid) 
structure and bubble formation. At 600 ps a bubble oc- 
cupies 25% of the confined space, and at 1 ns we find 
that 90% of the confined space is dewetted. With the 
size of the plate area that we simulate, it is not possible 
to observe simultaneously the vapor, liquid, and crystal 
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phases [only approximately five hexagons per side can be 
clearly observed in Fig. ^a)] . A precise investigation of 
the existence (or lack thereof) of a triple point requires 
larger plate surface areas than we have used in this work. 



B. Water structure: effect of pressure 

To study the effects of pressure on the structure of 
water we fix d = 1.6 nm. This separation is large enough 
so that water properties in the middle region between the 
plates are essentially those of bulk water ^| . 

Figure 13 shows the average water density between the 
plates at each pressure, p{P). For comparison, we show 
also the density of bulk water reported in [sl]- We de- 
fine p{P) = {nit)) X mH2o/Vconf where ■mH20 is the 
mass of a water molecule and {n(t)) is the average num- 
ber of molecules in the confined space between the walls. 
When dealing with confined systems, defining the ac- 
cessible volume for the confined liquid is not unique. 
We estimate Vconf in two ways: by the formal dimen- 
sions of the confined space, based on the defined d, i.e., 
Vconf = 1-6 X 3.217 X 3.215 nm^; and by an effective 
available volume Vconf = deff x3.217x3.215 nm^ where 
deff is an effective plate-plate distance. If Pn{z) is the 
probability density to find an O atom in a slab parallel 
to the plates located at z {J Pn{z) dz = 1), then deff 
is defined such that Pn{idef f /2) = 0.5 [see Fig. ISJ^a)]. 
The first approach to the calculation of Vconf underesti- 
mates p{P) because the formal confined volume extends 
practically to the plate surface (i.e., the plane where the 
H atoms would be located), and thus includes space not 
accessible in practice to water molecules. The second (ef- 
fective) definition of Vconf overestimates p{P) because it 
leaves out from the calculation a portion of volume that 
is in fact accessible to water molecules. 

p{P) in the confined space is always smaller than the 
bulk density for P < 0.06 GPa, independently of the 
method used to calculate Vconf- For P > 0.06 GPa the 
two estimates of Vconf produce drastically different re- 
sults. One expects that at high P the water density in 
the confined space approaches that of the bulk water. 
It can be seen from Fig. [7| that the formal definition of 
Vconf, using the d- value defined a priori, leads to a more 
reasonable high-P trend. For all cases, one finds that 
the slope of p{P) is always greater than that of bulk wa- 
ter indicating greater compressibility between the plates. 
To confirm this, we interpolate p{P) in Fig. with a sec- 
ond order polynomial, and calculate the corresponding 
derivative with respect to P. The resulting values for the 
compressibility, kt, at P = GPa and T = 300 K are: 
KT 0.47 GPa-\ for bulk water; and ht 1.47 GPa"! (for- 
mal volume) and kt 1-70 GPa~^ (effective volume), for 
confined water. The experimental value of the c ompress- 
ibility at these conditions is kt 0.45 GPa~^ |106| |. in good 
agreement with the simulations. 

Figure EI a) shows the density profile Psiab{z), i.e. the 
local density in a slab parallel to the plates, located at z. 



It was calculated using 39 slabs of width 0.0411 nm, with 
an area equal to the plates area. Psiab{z) resembles the 
number density at normal T and P reported in Ref. . 
Figure IHl^a) indicates that the main effect of increasing 
P is to push water molecules toward the plates. Accord- 
ingly, as P increases, Psiab{z) increases for z w and 
the maxima at z ± 0.51 nm for P = —0.05 GPa shift to 
z ± 0.63 nm at P = 0.2 GPa. Moreover, the maxima and 
minima of Psiab{z) are more pronounced at P = 0.2 GPa 
than at P = —0.05 GPa, as molecules increasingly sample 
'harder' plates at high P. We obtain the same conclu- 
sions when looking at the probability density function for 
the H atoms. 

To analyze the local packing and hydrogen bonding, we 
calculate the average coordination number CN{z) and 
the order parameter q{z) used to c haracterize the local 
tetrahedral order in water |l07l llOSj . Averages are taken 
over slabs of width 0.0411 nm centered at different z. We 
define CN{z) as the number of neighbor oxygens within 
a distance < 0.32 nm (the first minimum of the radial 
distr ibution function at p = 0.984 g/cw? and T = 28 4 K 
|109| ) from a central O atom. q{z) is defined as |107| 

334. ^.2 
9=1- g I] (cosVjfc + gj (4) 

where V'jfc is the angle formed by the lines joining the 
oxygen atom of a given molecule and those of nearest 
neighbors ] and fc (< 4). In this work, we also include 
the oxygen (O) atoms in the plates when considering the 
O nearest neighbors j and fc. 

FigurelHIb) shows CN{z) for different values of P. At 
all P, we find that CN{Z) for 2; « is close to 4 as is 
the case for bulk water. The slight increase of CN{z) 
with P at z = can be due to the fact that we are 
using a fixed cutoff of 0.32 nm to calculate CN (while the 
first minimum of the radial distribution function depends 
slightly on P). We observe that CN{Z) decreases as 
one approaches the plates. However, right next to the 
plates we find sharp maxima. Moreover, such maxima 
increase with P as water molecules are pushed toward 
the plates [see Fig. |Hta)]. When the O atoms in the 
plates are not considered in the calculation of CN{z), 
we find a distribution similar to that of In this case, 
CN{z) « 2.5 next to the plates, indicating that molecules 
in contact with the surfaces lose on average 1.5 — 2 water 
neighbors. 

Figure|SJc) shows q{z) for different values of P. At all 
P, we find that q{z) w 0.62 at z = 0, close to the value 
0.60 - 0.65 for bulk water at 0.85 < p < 1.15 g/cm^ and 
T = 300 K jl03 (indicated in the figure). The peaks of 
CN{z) next to the plates are accompanied by the rapid 
decrease of q{z) approaching the interface. At the plates, 
q{z) fa 0.5, and this value does not change with P, that 
is, the local tetrahedrality is not apparently affected by 
P for molecules in contact with the surfaces. 

To study the orientational structure, we follow and 
compute the distribution of angles 0hb between the four 
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HB vectors of water molecules and the inward pointing 
normal to the plates. Each molecule is associated with 
four HB vectors: these four vectors point tetrahedrally 
outwards from the oxygen Lennard-Jones site, and such 
that two of them join the O site to the same molecule's H 
atoms. The normalized distribution P{6hb) in the prox- 
imity of the plates is shown in Fig. |^for different values 
of P. P{9hb) shows maxima at Ohb — 70° and 180° indi- 
cating a preferred orientation for a water molecule with 
one HB vector pointing toward the nearest plate. These 
results, and a P{Ohb) calculated in slabs at different dis- 
tances from the plates (not shown), are in agreement with 
the results obtained in ref. in constant volume MD 
simulations. Interestingly, we find that when P increases, 
despite changes in Psiab{z), P{0hb) barely changes, as is 
evident in Fig. ^ In particular, molecules next to the 
plates preserve their distinct orientational order under 
pressure. This is in agreement with the noted invariance 
of q{z) « 0.5 next to the surfaces for all P. 

The radial distribution function (RDF) for the water 
molecule O atom and the O atom at the surface of the 
plate, go'o{r), is shown in Fig. ^| At all P, go'o[f) has 
maxima at 0.32, 0.60, and 0.77 nm. For comparison, the 
first maximum of the bulk water RDF at p = 0.984 g/cvo? 
and T = 284.5 K 109] is at 0.28 nm. Thus, water oxy- 
gen atoms are located farther from plate surface oxygen 
atoms than from other water molecule oxygen atoms, as 
expected for hydrophobic surface atoms. The increase of 
P does not shift the location of the maxima and minima 
of go'o{T), but changes their relative heights. As P in- 
creases, more molecules are found in the confined space 
and, at the same time, more molecules shift from inter- 
stitial to coordination shells. Thus, the confined liquid 
becomes more structured upon compression. 

In summary, the effect of pressure on the structure of 
water confined between hydrophobic plates is to enhance 
water structure, pushing the water molecules toward the 
plates. The orientation of water molecules next to the 
plates is not sensitive to compression, as indicated by 
Fig. m 



IV. RESULTS: HYDROPHILIC PLATES 

Over the range of conditions investigated here (T ~ 
300 K, P = -0.05, 0.05, and 0.2 GPa, and in each case 
d = 0.6, 1.0, and 1.6 nm), we find that the confined 
water remains in the liquid state, and shows no sign 
of capillary evaporation or crystallization. In agreement 
with the hydrophobic plate simulations, we find that for 
P < —0.1 GPa, the whole system cavitates; i.e., liquid 
water is unstable. Experimentally, this would manifest 
itself as loss of cohesion and the appearance of a macro- 
scopic vapor phase; computationally, it becomes impos- 
sible to maintain tensions in excess of 0.1 GPa and the 
system volume grows uncontrollably due to the appear- 
ance of the vapor phase. However, even in this case we 
find no capillary evaporation between hydrophilic plates. 



In other words, liquid water fills the confined space, and 
the hydrophilic plates induce wetting even when there is 
bulk cavitation around the plates. 



A. Water structure: effect of pressure 

As previously done for the case of the hydrophobic 
plates, we chose a distance between plates of d — 1.6 nm. 
The average density in the confined space between hy- 
drophilic plates is shown in Fig. ^2 As for the hydropho- 
bic plates case, we define p{P) = ("-(^)) x nT-Hio /Vconf 
and use both d = 1.6 nm and Pn{±deff) = 0.5 to com- 
pute Vconf- In contrast to the hydrophobic case, both 
methods now lead to very similar p{P) values. The un- 
derestimated (formal volume) and overestimated (effec- 
tive volume) values of p{P) are very close to each other 
and bracket those for bulk water. Interestingly, Fig. ITTI 
shows that the slopes of p{P) for confined water (using 
both volume definitions) are almost the same as that of 
bulk water. This suggests that the compressibility of wa- 
ter confined by hydrophilic plates is the same as that of 
bulk water. In fact, we calculate the compressibility as 
in Sec, nil 51 for the case of water confined by hydropho- 
bic plates. We find that, using either the formal volume 
or the effective volume, the compressibility of water con- 
fined by hydrophilic plates is k,t ~ 0.52 GPa^^, close to 
the value kt ~ 0.47 GPa^^ of bulk water. 

Figure El shows Psiab{z), CN{z), and q{z) for water 
confined between hydrophilic plates at different values of 
P. The three quantities arc insensitive to P over the 
range of conditions investigated here. Thus, the strong 
attraction exerted by the walls on liquid water is not 
measurably perturbed by the additional forces associated 
with compression. The two maxima of Psiab{z) in Fig. 
W^ Si) clearly indicate the presence of two water layers 
next to the plates. W e not e that Fig. ITW a') is very sim- 
ilar to Fig. 3 of Ref. TTol for the TIP4P water model. 
From Fig. (Elb) and (c), we find that both CN{z) and 
q{z) in the presence of hydrophilic plates decrease to zero 
very sharply at the plate surfaces. Moreover, a compar- 
ison of Fig. |S^c) and I12r c) shows that in the absence 
of the hydroxylated groups on the plates, the values of 
q{z) start to decrease at larger distances from the walls 
than in the case of hydrophilic plates. The resulting pro- 
files in Fig. I12r b) and (c) are almost flat for the region 
between the plates meaning that the values of CN{z) 
and q{z) next to the plates barely change with respect 
to their bulk water values at z = 0. These results sug- 
gest that water molecules far from the plates have simi- 
lar local environment (in terms of tetrahedral order and 
number of nearest-neighbors) as those next to the plates 
(when including the plate surface oxygen atoms). For 
the molecules next to the plates, the O — H atoms on the 
plate surface act as 'virtual' water molecules providing 
extra HBs to the real water molecules. Again, the slight 
change in CN{z) at z « is probably due to the fact that 
we use a fixed cutoff when calculating the coordination 
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number. 

FigurelTSlshows P{0hb) corresponding to a slab next to 
the plates of width 0.1 nm. P{9hb) is in agreement with 
the corresponding distribution computed in ref. |llCl| 
for the TIP4P water model and is complementary to the 
distribution found in the case of hydrophobic plates (see 
Fig. |5Jl . The presence of polarity on the plate allows wa- 
ter molecules to form HBs with atoms on the plate. The 
resulting P{9hb) increases abruptly for 9hb < 20° and has 
a sharp peak at 6hb ~ 110°. As already found in |llOj| for 
the TIP4P water model, these maxima can be explained 
by the presence of multiple HBs that a water molecule 
can have with Off groups associated with the plate (see 
Fig. 14 in ^3). Ultrafast electron crystallography mea- 
surements of water confined between hydrophilic walls 
[Tsf also show that water molecules interact at two sites 
of the substrate. 

The RDFs between the plate oxygen (O') or hydro- 
gen {H') atoms and the water molecule oxygen atoms 
are shown in Fig. ^1 These distributions do not change 
with P, confirming the insensitivity of confined water 
structure between hydrophilic plates to changes in exter- 
nal pressure. Figur e IT^ is in agreement with Figs. 11 
and 12 of Ref. [ll3, which were obtained for the TIP4P 
water model. 



V. SUMMARY 

We have presented results from MD simulations of 
water confined between nanoscale (hydrophobic or hy- 
drophilic) plates at T = 300 K and for a range of values 
of pressure P (-0.15 GPa < P < 0.2 GPa) and plate- 
plate separation d (0.4 nm < d < 1.6 nm). 

In the case of hydrophobic plates, a phase diagram 
in the P — d plane summarizing the MD results is pre- 
sented, and three phases (liquid, vapor, and bilayer ice) 
are identified. At low P we find capillary evaporation. 
The distance at which this drying occurs decreases with 
increasing P. Furthermore, the transition line in the 
P — d phase diagram separating the vapor and liquid 
phases is followed to small values of d, below which water 
molecules become individually trapped and immobilized 
by the surface atoms. It is possible that at higher T, this 
line ends in a critical point at high d (corresponding to 
the liquid- vapor critical point of bulk water). 

The bilayer ice is composed of two layers of hexagons 
in registry along the surface normal. The resulting ice 
is similar to that found in Ref. '^9l'|. This ice is found 
at P > 0.1 GPa and d « 0.6 nm, but only in a narrow 
range of d. The region corresponding to ice in the P ^ d 
phase diagram might be connected with the vapor-liquid 
transition line for larger plates, corresponding to a triple 
point. However, the plates we simulate are small, and 
we cannot observe simultaneously the liquid, vapor and 
crystal phases. Further investigation is required on this 



point. We also note that at T = 300 K and in the pres- 
ence of nanoscale plates, we do not find indications of a 
monolayer ice, like that observed in MD simulations us- 
ing^the TIP4P and TIP5P potentials and infinite walls 

When simulating water confined by hydrophilic (hy- 
droxylated silica) plates, we find that the confined wa- 
ter remains in the liquid phase at all P and d studied. 
Moreover, even at P where we observe cavitation in the 
bulk water, the confined water is in the liquid phase. In 
other words, the hydrophilic plates induce wetting in the 
confined space (a phenomenon reminiscent of capillary 
condensation). This result suggests that one way to sta- 
bilize liquid water under tension may be by hydrophilic 
confinement. Experiments of water confined between a 
hydrophobic and a hydrophilic surface (Janus interface) 
[73, suggest a similar physical picture, where the hy- 
drophobic surface encourages water to dewet, while the 
hydrophilic surface constrains water to be present. 

We also study the effect of P on the hydration of both 
hydrophobic and hydrophilic plates. We focus on the wa- 
ter molecule distribution, average coordination number 
and local tetrahedral order parameter along the direction 
normal to the plates. In the case of hydrophobic plates, 
all of these quantities indicate that as P increases, water 
molecules are pushed toward the plates while the liquid 
becomes more structured. However, the water molecule 
orientation next to the plates shows no dependence on 
P. As previously reported in ref. [4^ . molecules have on 
average one HB pointing to the plate. Strikingly, for the 
case of hydrophilic plates, we find no change in the liquid 
structure with P in the range —0.1 GPa to 0.2 GPa. 

The strong differences observed in the behavior of wa- 
ter confined between hydrophilic and hydrophobic plates 
leads to the interesting question of what the thermody- 
namic and structural properties of water are when it is 
confined between heterogeneous plates (e.g., hydrophobic 
plates that have been partially hydroxylated) . Results 
addressing this question will be reported in a subsequent 
report. 
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TABLE L Poten tial p arameters for plate-water interactions 
(taken from Ref. |lloj). 



Atom Type 


e [kJ/mol]("' 


a [nm]("' 


Charge [e]*'') 


O 


0.6487 


0.3154 


-0.71 


Si 


0.5336 


0.3795 


0.31 


H 






0.40 



Lennard- Jones parameter for plate-water O atom 
interactions (hydrophobic and hydrophilic plates). 

Charges on Si — O — H groups (hydrophilic plates). 



TABLE IL Simulation times for the different kinds of plates 
and plate separation d. 



Plate 


d [nm] 


Simulation time [ps] 


hydrophobic 


> 0.8 


200 


hydrophobic 


< 0.8^"' 


400'") 


hydrophilic 


< 1.6 


200 



For hydrophobic plates with d = 0.6 nm and 
P = 0.05 GPa, the simulation time is 1 ns. 



***** • 



* . • 



• • • 
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FIG. 1: (Color online) (a) Front and (b) lateral view of 
the fully hydroxylated silica plate [(1.1.1) octahedral face of 
cristobalite] , one of the two kinds of plates used in this study. 
White, red, and gray spheres represent hydrogen, oxygen, and 
silicon atoms, respectively, (c) y — z cross section showing the 
space available to the water molecules, (d) Top and (e) lat- 
eral view of the surface 5104 unit, in which the surface oxygen 
atom is hydroxylated. 
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FIG. 2: (Color online) Phase behavior of confined water as 
a function of pressure and separation between nanoscale hy- 
drophobic plates. The temperature is fixed at T = 300 K. 
Filled circles indicate liquid state points. Empty squares in- 
dicate vapor state points at which capillary evaporation oc- 
curs (i.e., we observe no water molecule between plates). For 
d « 0.6 nm, a bilayer crystal is formed as indicated by red 
X symbols. At d = 0.6 nm and P — 0.05 GPa, fluctuations 
between liquid and vapor states are observed for at least 1 ns. 
Continuous lines suggest a schematic phase diagram based 
on the simulation results. In the figure, the dashed red line 
corresponds to the assumption that the crystal region is not 
connected with the liquid- vapor transition line, i.e. that there 
is no 'triple point' between the liquid, the vapor, and the bi- 
layer crystal. The black dotted line is the extrapolation of the 
liquid- vapor transition line into the d < 0.4 nm region, where 
no liquid can be simulated due to the small distance between 
plates. 



FIG. 3: (Color online) Sequence of time-separated snapshots 
showing capillary evaporation between hydrophobic plates at 
P = GPa and d = 0.6 nm. Times correspond to (a) 18, (b) 
40, (c) 135, and (d) 239 ps. We only show water molecules in 
the slab corresponding to the space between plates (i.e., the 
plates themselves have been omitted for clarity) . The square 
cavity evident at 239 ps corresponds to the whole area of the 
plates. 
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FIG. 4: (Color online) (a) Front and (b) lateral snapshots 
of the bilayer ice formed between hydrophobic plates at P = 
0.2 GPa and d = 0.6 nm. 



d=0.5nm 
d=0.6nm 



d=0.4nm 



0.0-1 



100 200 300 

time [ps] 



400 




r [nm] 




UH ' 1 - ' : ' 1 - 

-0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 

z [nm] 



FIG. 5: (Color online) (a) Lateral mean square displacement 
(MSD) as a function of time at P = 0.2 GPa and different 
plate-plate separations d. Hydrophobic plates. At d = 0.6 nm 
the confined water crystallizes and the MSD for long times is 
smaller than for d = 0.5, 0.8 and 1.0 nm (where the system is 
in the liquid phase). For comparison, we also show the results 
for d — 0.4 nm where water molecules are trapped between 
the plate atoms due to the small value of d. (b) Radial distri- 
bution function (RDF) parallel to the plates, Qxyir), for the 
crystal (d = 0.6 nm) and liquid (d = 0.8 nm) phases. Inset: 
ratio of the crystal RDF to that of the liquid, (c) Probabil- 
ity density function associated with the distribution of water 
molecules between the plates (/ P„(z) dz — 1). 




FIG. 6: (Color online) (a) Schematic diagram showing the 
location of water molecules in the bilayer ice phase relative 
to the position of the hydrophobic plate surface atoms. Red 
filled circles represent 'oxygen' atoms belonging to a single 
plane on the plate surface. The next layers of 'oxygen' atoms, 
above and below this plane, are represented by empty and 
filled black diamonds, respectively. The atoms at the plate 
surface are arranged in tetrahedra pointing either into or out 
of the page of the figure [cf. Fig.0(d),(e)]. Tetrahedra point- 
ing into the page correspond to the cross-hatched triangles 
centered on the filled diamonds, while those pointing out of 
the page correspond to triangles centered on the empty dia- 
monds. Water oxygen atoms are represented by empty circles; 
a typical hexagon is indicated by blue lines. Lateral (b) and 
top (b) view of molecules belonging to the bilayer ice, showing 
the hydrogen bonds. We show only two water molecules with 
the corresponding hydrogen atoms (red and white spheres). 
These molecules belong to different layers of the ice and their 
corresponding nearest neighbors are represented by blue and 
green spheres. Yellow lines are a guide to the eye showing the 
lattice characterizing the bilayer ice. 
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FIG. 7: (Color online) Average density p{P) for different val- 
ues of P in the confined space between hydrophobic plates, 
and in bulk water. Circles correspond to the values of p{P) for 
bulk water taken from ref. |88^|. For confined water, we show 
two estimates of p{P) corresponding to different definitions 
of the confined volume. Squares and triangles correspond to 
p(P) calculations when using the 'effective' and 'formal' def- 
inition of volume, respectively (see Sect. Illl'Hl for details). 
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FIG. 8: (Color online) (a) Density profile, Psiab{z) (i.e., the 
local density in a slab of width 0.0411 nm, parallel to the 
plates, located at z); (b) average coordination number CN{z); 
and (c) order parameter q{z). Arrows in (c) indicate the range 
of q values for bulk water at T = 300 K and 0.85 g/cm^ < 
P ^ 1.15g/cm'^ 107j . Hydrophobic plates. 
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FIG. 9: (Color online) Normalized distribution P{9hb) of an- 
gles 8hb between the hydrogen bond vectors of water molecules 
and the inward pointing normal to the plates. Average is per- 
formed over molecules at a distance < 0.2 nm from the plates. 
Hydrophobic plates. 




FIG. 10: (Color online) Radial distribution function for the 
water molecule O and plate surface O atoms. Hydrophobic 
plates. 
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FIG. 11: (Color online) Density of water between hydrophilic 
plates compared to that of bulk water. Symbols and expla- 
nation are the same as in Fig.[7| 
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FIG. 12: (Color online) Same as Fig. |H1 but for the case of 
hydrophilic plates. 
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FIG. 13: (Color online) Normalized distribution P{6hb) of an- 
gles 9hb between the hydrogen bond vectors of water molecules 
and the normal direction to the plates, for molecules at a dis- 
tance < 0.1 nm from the plates (c/Fig.inj. Hydrophilic plates. 
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FIG. 14: (Color online) Radial distribution lunctions for the 
plate oxygen atoms (O'), and (a) the water oxygen and (b) 
water hydrogen atoms. (c)(d) Same as (a), (b) for the plate 
hydrogen atoms {H') . Hydrophilic plates. 



